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ABSTRACT 

We present 26 point-sources discovered with Chandra within 200" (« 20kpc) of the center of the 
barred supergiant galaxy NGC 1365. The majority of these sources are high-mass X-ray binaries, 
containing a neutron star or a black hole accreting from a luminous companion at a sub-Eddington 
rate. Using repeat Chandra and XMM-Newton as well as optical observations, we discuss in detail 
the natures of two highly-variable ultraluminous X-ray sources (ULXs): NGC 1365X1, one of the 
most luminous ULXs known since the ROSAT era, which is X-ray variable by a factor of 30, and 
NGC 1365X2, a newly discovered transient ULX, variable by a factor of >90. Their maximum X-ray 
luminosities (3-5x10^" ergs" ^, measured with Chandra) and multiwavelength properties suggest the 
presence of more exotic objects and accretion modes: accretion onto intermediate mass black holes 
(IMBHs) and beamed/super-Eddington accretion onto solar-mass compact remnants. We argue that 
these two sources have black-hole masses higher than those of the typical primaries found in X-ray 
binaries in our Galaxy (which have masses of < 20 Mq), with a likely black- hole mass of 40-60 in 
the case of NGC 1365 XI with a beamed/super-Eddington accretion mode, and a possible IMBH in 
the case of NGC 1365X2 with M = 80-500 M©. 

Subject headings: GALAXIES: INDIVIDUAL (NGC 1365), X-RAYS: binaries, X-rays: galaxies. 



I. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are bright com- 
pact off- n uclear point sources detected in nearby galax- 
ies fe.g.. iMakishima et"an 120001: iFabbiano et"al] l2001l : 
iSwartz et al.ll2004fl . For sub-Eddington accretion rates in 
the absence of beaming, their bolometric luminosities of 
a fewx 10'^^-10'*° erg s^^ require black hole masses A/bh ^ 
10-400 M0, implying the existence of intermediate-mass 
black holes (IMBHs). Despite increasing solid evidence 
for the existence of supermassive bl ack holes in sizable 
bulge galaxies (Mbh > 10^ Mq; e.g.. iGreene fc Holl2007l . 
and references therein) and the numerous examples of 
stellar-mass black holes (SMBHs) in our Galaxy and 
beyond, there are few convincing candidates for black 
holes with masses in the intermediate regime. The term 
SMBHs, as used here, refers to black holes produced 
at the end stage of normal stellar evolution. Obser- 
vationally, SMBHs are constrained to have masse s be- 
low - 20 M© fe.g. iRemiUard fc McChntockl l2006l and 
references therein; see also Fryer fc Kalogera 2001 and 
Orosz et al. 2007). Theoretical studies show that com- 
pact remnants of up to 60-80 M© could form under spe- 
cial circumstances (assuming solar abunda nces and ini- 
tial masses of up to about 200 Mq; e.g., iHeger et al.l 
120031: lYungelson et al.l I2008D and we adopt this lower 
limit^ as our definition of an IMBH: a compact rem- 
nant with mass between 100 and 10^ M©. The existence 
of IMBHs presents challenges in terms of formation and 
feeding mechanisms leading some researches t o prefer al- 
ternative explanations for ULX emission fe.g.. lKhIill2006l : 
lYungelson et al.ll2008[ ). 

^ Max Planck Institute for Extraterrestrial Physics, Postfach 
1312, Garching, 85741, Germany 

^ For metallicities much below solar, characteristic of the early 
universe, black holes with m asses over a few hundred solar masses 
can form HHeger et al.|[200l) . 



A simple explanation for the observed high X-ray lu- 
minosity of an ULX is a misidentified background ac- 
tive galactic nucleus (AGN), whose chance coincidence 
with a nearby galaxy and low optical emission contrive 
to suggest a more unusual object. A recent example 
is one of the most lumin ous ULXs found to date in 
the galaxy MCG-03-34-63 (jMiniutti et al.ll2006l ). which 
was was later identified with a background AGN (Mini- 
utti 2008, private communication). A more intriguing 
alternative for the origin of the ULX emission than a 
background AGN is super-Eddington accretion from a 
high mass X-ray binary (HMXB) with a SMBH pr i mary . 
If the emission is beamed (e.g.. [ Pou tanen et al.l 120071 : 
iKording et al.ll200l : lGeorganopoulos et al. 2003) and the 
outflow or jet axis points close to out line of sight, emis- 
sion of up to ~ lO'^'^ergs^^ can be explained without 
the need to invoke IMBHs. Super-Eddington and su- 
perluminal accretors have been observ ed in our Galaxy 
(e.g. y4641, SS433, GRS 1915 + 105: iRevnivtsev et alJ 
120021: iBegelman etl l 2006; Mir abel fc Rodriguezlll99l 
and according to [King (2006) could be responsible for 
the apparent super-Eddington emission from the major- 
ity of ULXs. Sources with X-ray luminosities in excess 
of 10'*'^ ergs~^, however, are much rare r and m ore likely 
to host IMBHs. Indeed ISwartz et all ([2004') find that 
the ULX luminosity function for star-forming galaxies 
requires a broken power-law fit with break luminosity of 
10^°ergs~^, consistent with the idea of a change in the 
nature of the ULX population above this luminosity. 

Luminous ULXs are more commonly found in star- 
forming galaxies, leading researchers to suppose that the 
ULX phenomenon is associated with recent star forma- 
tion and p ointing towards a HMXB origin for the ULX 
emission. iRappaport et al.l (|2005D present a detailed 
analysis of HMXB evolution models together with a bi- 
nary population synthesis code to study the theoretical 
expectations for ULX populations of star-forming galax- 
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TABLE 1 
NGC 1365: Chandra X-RAY SOURCES 



Source 


RA 


Dec 


Fo.3-10kcV 


r 


Norm 


C-stat 










(.5; 




(') 


X-3 


03:33:21.15 


-36:08:15.1 


1 g+3.6 


1.7±0.7 


2.73E-6 


18/13 


X-4 


03:33:23.11 


-36:07:53.1 


9 O+2.0 


1.8±0.3 


3.71E-6 


26/25 


X-5 


03:33:26.03 


-36:08:38.3 


9 9+3.2 
^■^-1.8 


1.5±0.5 


2.57E-6 


17/20 


X-6 


03:33:26.28 


-36:08:49.6 


9 9 + 2.2 
^■^-1.9 


1.4±0.4 


2.20E-6 


18/18 


X-7 


03:33:26.64 


-36:08:13.9 


1 5+1-'' 


2 fixed 


2.61E-6 


2/7 


X-8 


03:33:29.63 


-36:08:29.9 


^■0-1.2 


1.9±0.5 


2.77E-6 


17/19 


X-9 


03:33:29.64 


-36:07:56.1 


U-8_0.8 


2 fixed 


1.48E-6 


11/8 


XIO 


03:33:30.38 


-36:08:21.5 


U-8-0.8 


2 fixed 


1.48E-6 


11/7 


Xll 


03:33:31.09 


-36:08:07.7 




O.litO.9 


6.76E-7 


17/10 


X12 


03:33:31.59 


-36:08:08.7 




1.4±0.3 


5.82E-6 


34/46 


X13 


03:33:32.22 


-36:06:42.8 


4.61^:^ 


1.1±0.5 


3.26E-6 


16/22 


X14 


03:33:32.37 


-36:09:02.9 


3 0+'^° 
•^■"-3.0 


0.8±0.8 


1.37E-6 


9/10 


X15 


03:33:34.16 


-36:11:02.5 




1.5±0.2 


1.21E-5 


61/83 


XI 


03:33:34.61 


-36:09:36.6 


i4t^:i 


l.SitO.l 


1.69E-5 


90/104 


X16 


03:33:36.45 


-36:09:58.0 


^■'^-1.0 


1.1±0.5 


7.32E-7 


6/5 


X17 


03:33:37.25 


-36:10:26.2 


9 q+7.0 


1.1±0.7 


2.14E-6 


10/10 


X18 


03:33:38.02 


-36:09:35.2 


S.9tli 


1.4±0.2 


l.OOE-5 


60/75 


X19 


03:33:39.21 


-36:10:01.9 


-1 q + 2.0 

^"^-1.3 


1.5±0.4 


1.66E-6 


14/13 


X20 


03:33:39.75 


-36:10:37.3 


i.stli 


2.4±0.5 


3.80E-6 


17/21 


X21 


03:33:40.28 


-36:07:26.7 


4.5t|l 


1.5±0.4 


5.11E-6 


44/32 


X22 


03:33:41.94 


-36:07:41.3 


6 0+^° 
"■^-6.0 


-0.1±0.5 


6.32E-7 


8/14 


X23 


03:33:44.48 


-36:09:58.3 




1.2±0.5 


3.40E-6 


15/15 


X24 


03:33:49.65 


-36:10:25.6 


9 7+1.7 
^- ' -1.4 


1.8±0.3 


4.53E-6 


23/31 


X25 


03:33:37.33 


-36:07:23.1 




1.1±0.4 


1.75E-2 


22/20 


X2 


03:33:41.85 


-36:07:31.4 




0.5±0.3 


2.47E-6 


47/44 


X26 


03:33:44.29 


-36:07:17.7 


rj 7 + 3.2 

'■'-2.8 


1.2±0.2 


6.48E-6 


75/68 



Note. — All but the three sources below the line were detected in our 2002 
Chandra observation (ObsID 3554); the last three sources are from the 2006 
Chandra observations with ObsID 6871. Column (1) gives the source designa- 
tion; (2) and (3) the source position, J2000; (4) the Galactic absorption cor- 
rected 0.3-10 keV flux in units of 10^^** ergs~^ cm~^; (5) the power-law index 
of the C-stat fit; (6) the 1 keV normalization of the C-stat fit in units of pho- 
tons keV~^ cm~^ ; (7) the C-stat per degree of freedom. 



ies under the assumption of a HMXB origin. They con- 
clude that the majority of ULXs with ^ 10^°ergs~^ 
can be explained in terms of binary systems containing 
SMBHs. Consequently, from both theoretical and empir- 
ical perspective, sources with X-ray luminosities in ex- 
cess of lO'^'^ergs"^ remain the most puzzling examples 
of ULX emission. 

In this paper we characterize the X-ray point-source 
population of the giant barred spiral galaxy NGC 1365. 
We show that it is compatible with the population of high 
mass X-ray binaries expected for gal axies with the star - 
formation rate of NGC 1365. K omossa fc Schulj ()1998f ) 
discussed the first ULX in NGC 1365, NGC 1365X1, us- 
ing ROSAT HRI dat a. The source was also noticed by 
llvomoto et al] ()1997D in an Advanced Satellite for Cos- 
mology and Astrophysics (ASCA) observation. At the 
time of discovery this was one of the most luminous ULXs 
known, which, together with the unusual properties of 
NGC 1365 itself, motivated our subsequent 2002 Chan- 
dra follow-up observation. In this paper, we focus our 
attention on the X-ray flux and spectral variability of 
the two brightest ULXs in NGC 1365, which reach X- 
ray luminosities in excess of 10''°ergs~^, and are among 
the most luminous ULXs known to date (for compari- 



son samples, see iMakishima et a l."2000': ' Fabbiano et al.l 
l200lHSwartz et al.ll2004l; iDevi et al, 20071 Adding opti- 
cal and infrared (IR) data to the detailed X-ray modeling 
of Chandra and XMM-Newton observations, we evalu- 
ate the plausibility of the different emission scenarios 
for the two ULXs. The paper is organized as follows: 
In § [21 we present the X-ray point-source population 
of NGC 1365. The X-ray spectral and flux variability 
of NGC 1365X1, together with the available optical/IR 
data are analyzed in § [3l We discuss the X-ray and opti- 
cal data of the new ULX, NGC 1365X2, in §|4l and study 
in detail the plausibility of the different emission mecha- 
nisms for both ULXs in § [5l followed by a short summary 
in § [6l Throughout this paper we assume a luminosity 
distance of Dl = 21.2 Mpc to NGC 1365, and the stan- 
dard concordance cosmology of S pergel et al.l (|2003( ) . All 
Chandra and XMM-Newton data used in this paper were 
reprocessed using CIAO v.3.4 and HEASOFT v.6.2 (SAS 
v.7.1) starting with the event 1 files for Chandra and re- 
peating the basic pipeline processing for XMM-Newton 
followed by the standard calibration and filtering oper- 
ations in each case. All quoted flux and luminosity un- 
certainties are 90% confidence limits, unless noted oth- 
erwise. 
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Fig. 1. — Left: The background-subtracted, Gaussian-kernel-convolved image of the inner 6' X 6' 36kpc X 36kpc) of NGC 1365 in 
the 0.3— lOkeV band during our 2002 Chandra observation. North is up and East is to the left. The 26 X-ray point sources, detected in 
Chandra observations are denoted by 6"-radius circles. The two ULXs studied in detail here are marked by XI and X2. The two point 
sources marked "ON" as well as X2 were much fainter and not detected during the 2002 exposure shown here; the point source marked 
"OFF" fell below the detection limit between 2002 and 2006. Right: The cumulative number of X-ray point sources with 0.3-10 keV 
luminosity greater than a given value for the 23 sources detected during the 20 02 exposure (open histogram) and the expect ed number of 
background sources based on the logA'^ — log 5 results of lHasinger et aLl 11200 II) (hatched histogram). The dashed line is the I Grimm et aLl 
II2006I ) prediction for a spiral galaxy with SFR=12M,7- ^"--^ 
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2. THE X-RAY POINT-SOURCE POPULATION OF 
NGC 1365 

NGC 1365 (RA, Dec 03:33:36.40, -36:08:25.7; dis- 
tance modulus'^ of 31.6 mag) is a luminous giant barred 
spiral galaxy with a high star form ation rate (SFR ~ 
12M ( T,yr-^ for = 6.8 x lO" L^: [Lonsdale fc Heloiil 
119851 iKennicutt) Il998f ). Motivated by the discovery of 
NGC 1365X1 with ROSAT, we observed NGC 1365 with 
Chandra in 2002, in order to follow up the luminous ULX 
and to characterize the X-ray point-source population of 
the galaxy. The 2002 Chandra observation shows dif- 
fuse extended nuclear emission and 23 off-nuclear point 
sources'* in the central ^ 6' x 6' region. The point sources 
were detected usin g the wavdetect sl^OT ithm with source 
significance > 3(7 (iFreeman et al.ll2002f) as implemented 
in CIAO v.3.4. Spectra were extracted for each off- 
nuclear source from a 2"-radius aperture and fit with 
a simple power-law model assuming no intrinsic absorp- 
tion above the Galactic value using Cash-statistics. The 
positions, spectroscopic-fit parameters, and model-flux 
estimates of all off-nuclear point sources are presented 
in Table [Ij Between 2004 and 2006, NGC 1365 was the 
target of extensive XMM-Newton (P.I. G. Fabbiano) and 
Chandra (P.I. G. Risaliti) campaigns aimed at unravel- 
ling the highly obscured and variable active galactic nu- 
cleus with Chandra and studying the ULX NGC 1365X1 
with XMM-Newton. 

The left panel of Figure [1] shows the 2002 Chandra 
background-subtracted and smoothed image of the in- 
ner 6' X 6' (~ 36kpc X 36kpc; compare this to the 
B = 25™ arcsec^^ isophotal radius of 3Gkpc). The cu- 
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Database, 



|http: / /nedwww. ipac.caltech.edu/ 

* Three additional point sources were detected in later Chandra 
observations of the same area; all three are strongly variable and 
were too faint to be detected during the 2002 Chandra observation. 



mulative number of X-ray point sources with luminosities 
greater than a given value are shown as a histogram in 
the right panel. The X-ray luminosities of local star- 
formi ng galaxies ar e proportional to their SFRs (see 
,Gilfanov et al.|[2004 and references theirin) as the num- 
bers of HMXBs, which dominate the extended emission 
in local star-forming galaxies, increases wit h increasing 
SFR. Correcting for the SFR, I Grimm et all (|2006) have 
derived a simple analytical expression which predicts the 
differential number of X-ray point sources N for a given 
point-source luminosity, L38 (the luminosity in units of 
1038ergs-i): dN/dLss = 3.3 x SFR x Lg/'^^ We plot 
the integral form of this relation in the right panel of Fig- 
ure [T] with a dashed line. The expected number of X-ray 
point sources in NGC 1365 agrees well with the observed 
histogram for bright sources (logLx > 38.5), suggesting 
that the majority of the X-ray point sources are likely 
associated with the galaxy. The shaded histogram in the 
right panel of Figure [T] is a rough estimate of the possible 
contamination from background AGNs for each luminos- 
ity bin, which we compute using the minimum 2-10 keV 
luminosity for e ach b in and the logN — log S relation of 
Hasi nger et"aLl (|2001l ) as detailed in § [5TT1 We expect 
at most 2 of the 16 sources brighter than lO'^^ergs"^ in 
the 6' X 6' area to be background AGNs. We mention 
in passing that there is no correlation between the mea- 
sured point-source fluxes and their projected distances 
from the center of NGC 1365. 

The positions of the 26 X-ray point sources on an op- 
tical image of NGC 1365 are shown in Figure [S] Despite 
the fact that many of the X-ray point sources coincide 
with the minor and major spiral arms, only 4 out of the 
26 have an optical counterpart within ^ 2". Three of 
the four optical counterparts are extended sources, most 
likely H II regions, while the remaining southernmost 
source with an optical counterpart (XI 5 in Table [T|) is a 
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Fig. 2.— VLT FORS2 image of the inner 6' X 6' (~ 36kpc X 36kpc) of NGC 1365. 
denoted by 3"-radius circles here. North is up and East is to the left. 



The 26 X-ray point sources shown in Figure [Tl are 



point source. The optical counterpart (in the case of X2) 
and upper limits (in the case of XI) are discussed in § [j] 
and § 13. 3[ respectively. The remaining three sources with 
optical counterparts, X17, X15, and X13 (see Table [1]) 
have X-ray fluxes ofFo.3_iokov — 2.9x lO^^^'ergs"^ cm^^, 
-Fo.3-iokcv = 9.8 X 10""ergs"^ cm"2, and Fo.a-iokcv = 
4.6 X lO^^'^ergs"^ cm~^, and power-law spectral indices 
of r = f.l±0.6, r = 1.5 ±0.2, and T = 1.1 ±0.5 (assum- 
ing absor ption equal to t he Galactic value). The USNO- 
Bl.O (|Monet et al.ll2003D counterparts of X17 and X15 
have R-band magnitudes of Rl = 18.4 and R2 = 19.2 
(X17) and Rl = 19.2 (X15). 

3. NGC1365X1 

NGC 1365X1 (03:33:34.61, -36:09:36.6) is the most 
luminous ULX in NGG 1365. It was discovered by 
iKomossa fc Schuli (|1998l ) to be highly variable on long 
timescales (a factor of 10 within 6 months), with a peak 
luminosity reaching 5 x 10'"' erg s~^ in the 2-10 keV 
band (an ASCA measurement by llvomoto et al.l [T997t 
updated to Dl = 21.2 Mpc; note that ISoria et all (|2007l ) 
find a somewhat lower value by reanalyzing the ASCA 
spectrum). XI is ^ 74" south from the NGC 1365 nu- 



TABLE 2 
NGC1365: Observation Summary 



ObsID 


MJD 


XI Counts 


X2 Counts 


Tcff 


(1) 


(2) 


(3) 


(4) 


(5) 


Chandra 




ACIS-S 


ACIS-S 


ACIS-S 


3554 


52632 


167 


< 4 


9.9 


6871 


53835 


613 


51 


13.4 


6872 


53837 


880 


364 


14.6 


6873 


53840 


886 


306 


14.6 


6868 


53843 


745 


208 


14.6 


6869 


53845 


541 


117 


15.5 


6870 


53848 


190 


535 


14.6 


XMM 


pn/MOSl/MOS2 




pn/MOS 


0151370101 


52655 


151/70/84 




14.4/17.6 


0151370201 


52679 


24/18/12 




2.5/5.4 


0151370701 


52864 


70/15/30 




6.6/8.2 


0205590301 


53021 


1519/490/534 




49.0/56.5 


0205590401 


53210 


999/505/461 




30.0/44.6 



Note. — (1) Observation ID; (2) Modified Julian Date 
(MJD); (3) XI counts in the 0.3-10 keV band; (4) X2 counts 
in the 0.3— lOkeV band; (5) The eff'ective exposure time in units 
of 1000 seconds, after corrections for flaring and chip position. 
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Fig. 3.— NGC1365X1: light curves in the soft (0.2-5 keV; top 
panels) and hard band (5—10 keV; middle panels) and the fraction 
of hard (5— lOkeV) vs. total (0.2— lOkeV) emission as a function of 
time (bottom panels). The solid circles denote the ACIS-S Chan- 
dra count rates (CRs) while the open circles denote the ACIS-S 
CRs equivalent to the observed XMM-Newton EPIC-pn CRs. The 
error-bar propagation assumes Gaussian distribution of the uncer- 
tainties and Poisson errors in the measured counts. The right pan- 
els show the Chandra flare of April 2006. 
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Fig. 4. — NGC1365X1: luminosity variation with time. The 
triangles show the ROSAT / ASCA data, with error bars indicating 
only the uncertainty due to the assumed power-law model. The 
Chandra (XMM-Newton) power-law model estimated luminosities 
are shown with solid (open) circles with errorbars indicating the 
model-fit uncertainties; the inset zooms on the April 2006 Chandra 
flare. 



cleus (equivalent to a projected distance of w 7.5 kpc) 
behind an average Galactic column of 1.4 x 10^° cm"^. 
It is one of the brightest ULXs known, and one of the 
few ULXs with extensive X-ray (7 Chandra and 5 XMM- 
Newton pointings since 2002, in addition to older ASCA 
and ROSAT data) as well as extensive optical/IR cov- 
erage from the Very Large Telescope (VLT), the Hubble 
Space Telescope (HST), and the Two Micron All Sky 
Survey (2MASS). NGC1365 was the target of a recent 
Chandra campaign to observe the rapid changes of its 
obscured active nucleus. Table [2] below gives a summary 
of these 6 Chandra pointings as well as the 2002 Chan- 
dra and 2003-2004 XMM-Newton (5 pointings) observa- 
tions. We reprocessed all data using CIAO v. 3. 4. 1.1 and 
HEASOFT V. 6.1.2; the effective exposure times of each 
pointing, after corrections for flaring and chip positions, 
are listed in the last column of Table [D 

3.1. XI: X-ray Variability 

The soft (0.2-5 keV) and hard (5-10 keV) band hght 
curves of XI are shown in Figure [31 with the right panels 
focusing on the recent Chandra data. The ULX source 
underwent a rapid brightening in both soft and hard 
bands with the soft-band count rate apparently slightly 
leading the hard band. After varying by less than a fac- 
tor of ~3 (5) in the soft (hard) band over 2 years, the 
ULX source flared over the 13-day Chandra observation 
period, changing by a factor of ~5 (10) in the soft (hard) 
band. Most of the observed counts in both the Chandra 
and XMM-Newton exposures are in the soft band below 
5keV. The bottom panels of Figure [3] show the variation 
of the fraction of hard band photons (using equivalent 
ACIS-S counts) with time. Note that the source ap- 
pears to harden during the Chandra flare, but as a whole 
the contribution of the 5-10 keV emission to the total is 
small, typically 2-4%. Despite XMM-Newton' s higher ef- 
fective area above 5 keV a factor of 3.5 more counts for 
an EPIC-pn MED filter at the same flux level) there are 
too few counts in all but the two longest XMM-Newton 
observations above 5keV to obtain good constraints on 
any hard-band components. We return to this point in 
§[S21 below. 

Figure m shows the long term 0.5-10 keV X-ray vari- 
ability of NGC 1365X1, including the earlier ROSAT 
and ASCA observations.^ The X-ray luminosity of 
NGC 1365X1 is typically a few times lO'^^ ergss"-^, but 
varies between 2 x 10^^ ergs s~^ and 4 — 6 x 10*° ergs s~^ , 
or a factor of < 30 over an eight year period. Note that 
this large variation on a timescale of years is coupled 
with much quicker (on the order of days) variation by 
a factor of 6 during the April 2006 Chandra flare. We 
searched the longest (with non flaring intervals > 10 ks) 
XMM-Newton and Chandra exposures for variability on 
scales of minutes to hours but found no statistically sig- 
niflcant changes (the probabilities of constant rates were 
larger than 10% in all cases). This is consistent with pre- 



We 



iKo 



used t h e RO SAT HRI count rates given in 
; Schuld (HUD and the average power-law model 
fit obtained in § [32] below, F = 1.8, TYh,! = 10^1 cm-^, which 
agrees reasonably well with the earlier ASCA and PSPC model 
fits. The errorbars for the ASCA and ROSAT points indicate 
only the uncertainties associated with adopting the specific power- 
law model, ignoring measurement errors and the uncertainties 
associated with the use of different instruments and passbands. 
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Fig. 5. — NGC 1365 XI: variability of the spectral parameters for the absorbed power-law fits with MJD (left panels) and 2-10 kcV band 
flux (right panels). 



vious results on samples of ULXs, the majority of which 
show no short term variability (e.g., only 5- 15% of Chan- 
dra U LXs are variable on short timescales; iSwartz et al.l 
[200l . 

3.2. XI: Spectroscopic Fits 

We extracted spectra from 2"-radius (Chandra) or 15"- 
radius (XMM-Newton) regions'' around NGC 1365X1 
and selected background regions free of point sources at 
similar distance from the NGC 1365 nucleus to minimize 
contamination. The spectral fits were performed using 
Cash-statistics for individual spectra with less than 300 
counts (two Chandra and three XMM-Newton observa- 
tions marked "C-stat" in Table [3]) using 15-count binning 
for the remaining five Chandra and 20-count binning for 
the two longest XMM-Newton observations. In order 
to test specific models of X-ray emission we attempted 
to fit all spectra with two simple models: (1) power- 
law fits, including intrinsic absorption; (2) multi-color 
disk (MCD) plus power-law emission models with ab- 
sorption equal to the Galactic value. The second model 
is physically motivated, representing the optically thick 
thermal X-ray emission from the accretion disk {diskbb 
model in XSPEC) as well as a harder, corona comp- 
tonizcd component (represented by a simple powerlaw 
model in XSPEC). The spectral fits for each of the Chan- 
dra and XMM-Newton observations for the two different 
models are given in Tables |3] and HI respectively. 

The addition of an intrinsic absorber substantially im- 
proved the power-law model fits for 9 of the 12 observa- 
tions shown in Table [3l For the remaining three observa- 
tions with limited photon statistics, we fix the intrinsic 
absorption to the value estimated for the well constrained 

^ The size of the Chandra extraction region corresponds to over 
90% encircled energy. For XMM-Newton the 15"-radius extraction 
region was selected to minimize contamination from the obscured 
active nucleus in the case of XI and the encircled energy is expected 
to be about 65-70%. We extracted spectra from larger radii (up to 
45") for the two longest XMM-Newton exposures and found that 
the flux estimates using these larger radii are consistent with the 
ones quoted, within the errorbars. 



observation closest in time. The variability of the spec- 
tral parameters for the power-law fits are shown in Fig- 
ure [5l The required intrinsic absorption, iVn^i, ranges be- 
tween (7±2) and (21±3) x 10^" cm~^, a variation by a fac- 
tor of < 3. Similar increase in the absorbing column with 
increasing fiux wa s observed in the HMXB C ygnus-X3 
with INTEGRAL (Hialmar sdotter et al.l[200l . The es- 
timated power-law photon indices vary between 1.6 ±0.1 
and 2.5 ± 0.2, with the majority of fits having F = 1.9-2. 
From Figure [5l the spectral index appears to flatten with 
increasing flux during the Chandra flare of April 2006 
(indicating a harder spectrum) while the intrinsic ab- 
sorption appears to remain constant over the same 13- 
day period. Note that the spectral parameter variation 
in general needs to be interpreted with caution, due to 
the known degeneracy between F and Nm in power- 
law fits, as well as the lower sensitivity of Chandra at 
high energies (see also the d iscussion in Appendix A3 of 
'Portcgics Zwar t et al.ll2004[ ). 

Figures [6] and [7] show the Chandra and XMM-Newton 
fits for the second model, representing the direct accre- 
tion disk emission plus corona-comptonized component. 
Statistically these fits are indistinguishable from the ab- 
sorbed power-law fits presented above (compare the last 
columns of Tables[3]and[4|) . In fact if we only consider the 
0.3-5 keV emission, the MCD fits are adequate and no 
power-law component is necessary in any of the 12 obser- 
vations. Considering the full 0.3-10 keV band, the MCD 
model alone provides good fits to all but the three longest 
XMM-Newton observations, for which the addition of a 
power-law component improves the fit with >99% confi- 
dence. As noted in ij 13. li the Chandra and (to a lesser ex- 
tent) XMM-Newton effective areas result in insufficient 
counts above 5keV for all but the two longest XMM- 
Newton exposures. The additional power-law component 
required by these spectra has a photon index of F « 1. 
Fitting the longest exposure, XMM-Newton 0205590301, 
we obtain the normalization of this additional power law 
component and proceed to include a F = 1 power-law 
component in all MCD fits, scaling the normalization of 
the power law according to the observed 5-10 keV count 
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TABLE 3 

NGC 1365X1: Absorbed Power-Law Fits 



ObsID 


.^2-10 kcV 


.^0.5-2 koV 


r 


Nh 


Data/Method 


stat/DoF 


('-) 


\/^) 








I"; 


{' ) 


3554 




7 4+2.2 
'■^-4.9 


2.2±0.2 


20±7 


ACIS-S/C-stat 


80/103 


6871 


16.5t?| 




1.9±0.1 


11±4 


ACIS-S/x^ 


31/34 


6872 


34.0+^-^ 


24.61^:? 


1.9±0.1 


20±3 


ACIS-S/x^ 


36/49 


6873 


42.7lt-g 


19.0+2g 


l.fiitO.l 


13±3 


ACIS-S/x^ 


50/48 


6868 


35.4l^J 


l'-^-4.6 


1.6±0.1 


14±4 


ACIS-S/x^ 


36/42 


6869 


20.11^8 


12 4+2-'^ 


1.8 ±0.2 


13±4 


ACIS-S/x^ 


18/30 


6870 


5 1+1-5 
0-i_1.4 


"■1-1.6 


2.0 ±0.1 


13 fixed 


ACIS-S/C-stat 


101/120 


0151370101 


^•"-2.0 


1 7+0.6 
1- ' -0.7 


1.6±0.2 


8±5 


PN+MOS/C-stat 


300/332 


0151370201 


1 5+0.8 
'-■^-0.7 


3.2+^-7 


1.6±0.2 


8 fixed 


PN+MOS/C-stat 


89/74 


0151370701 


^•"-0.6 


1 c+o:9 

1-^-0.6 


2.5±0.2 


10 fixed 


PN+MOS/C-stat 


151/183 


0205590301 


0-^-0.8 


"■1-0.6 


2.06±0.08 


10±2 


PN+MOS/x^ 


144/142 


0205590401 


7 1+1-1 

'•1-1.1 


fi 7+1-° 
"■ '-0.8 


2.06±0.08 


7±2 


PN+MOS/x^ 


107/107 



Note. — The fluxes are quoted in units of 10"^'* ergs~i cm"^. All flux errors are 
90% confldence limits; all other parameter errors are la errors. (1) Observation ID; (2) 
Absorption corrected flux in the 2— lOkeV band; (3) Absorption corrected flux in the 0.5— 
2 keV band; (4) Power-law photon index; (5) Intrinsic absorption column density in units of 
lO^o cm-2; (6) Data and fitting method used; (7) x^/DoF or C-stat/DoF of the fit. 



TABLE 4 
NGC 1365X1: MCD-Power-Law Fits 



ObsID 
(1) 


.^diskbb 

(2) 


-Po.5-2 keV 

(3) 


.^2-10 kcV 

(4) 


(5) 


K 
(6) 


N 
(7) 


stat/DoF 
(8) 


3554 


4.6+?-3 


4 9+0.2 


g 0+0.3 


0.63±0.15 


0.014±0.011 


6.5E-6 


83/106 


6871 


^"-2.4 


19+0.5 
1^-2.5 


19+1-2 


0.86±0.07 


0.017±0.005 


8.9E-6 


31/35 


6872 


25 — 2.2 


i4in 


+ 11 


1.13±0.10 


0.0072±0.026 


17E-6 


49/50 


6873 


90 + 1.0 
^•^-4.0 




41+0-9 
^1-7.0 


1.17±0.13 


0.0054±0.0027 


23E-6 


49/49 


6868 


99+0.9 
^^-6.2 


i2lJ:^ 


34t^i 


1.22±0.13 


0.0046±0.023 


17E-6 


34/43 


6869 


,4+0.9 
1^-3.0 

6.5+0-? 




20«-I 


0.97±0.11 


0.0072±0.0033 


llE-6 


18/31 


6870 


3.4i°:| 


5.6+0-? 


0.95±0.13 


0.0036±0.0025 


2.2E-6 


105/120 


0151370101 


9 (,+0.3 
^•"-1.0 


i.8i;-? 


4.6t?-? 


1.02±0.20 


0.0012±0.0014 


2.7E-6 


297/333 


0151370201 


1 a+0.3 
l-«-0.7 


-1 9+0.3 
1-^-0.8 


0+0.3 

0.0_Q 7 


1.04±0.52 


0.0007±0.0023 


1.9E-6 


88/74 


0151370701 
0205590301 
0205590401 


2.8«-.g 
6.0^1-° 

6.7t^;^6 


2 0+^-2 

4 5+0-2 

c: + 0.2 
O-J_0.4 


2 0+°-2 
^•"-0.3 
7 + 1.9 
'•■^-2.1 
S 1 +0-2 
8-1-0.3 


0.46±0.07 
0.65±0.04 
0.52±0.03 


0.028±0.017 
0.015±0.003 
0.041±0.008 


1.3E-6 
4.4E-6 
5.5E-6 


151/183 
151/142 
99/108 



Note. — The fluxes are quoted in units of 10 l-^ergs 1 cm 2. (1) Observation ID; (2) Flux 
of the MCD component in the 0.0001-10 keV band; (3) Total flux in the 2-10 keV band; (4) Total 
flux in the 0.5— 2keV band; (5) Maximum MCD-flt temperature in keV; (6) Normalization of 
the MCD component, K = {Ri„/D)^ cos9, where i?i„ is the apparent inner accretion-disk radius 
in km, D is the distanc e to the source in units of lOkpc and is the inclination of the disk; 
see IKubota et aP l|1998l ): (7) IkeV normalization of the F = 1 power-law component in units 
of photons keV~l cm~2 g-l; (g) y^'^/DoF or C-stat/DoF of the fit; the type of the fit for each 
observation (x^ or C-stat) is the same as that given in Table[3] We assume no intrinsic absorption 
above the Galactic value of 1.4 X 1020 cm~2 for all fits. All flux errors are 90% confldence limits; 
all other parameter errors are la errors. 



rate.^ The resulting MCD plus power-law fit parame- 
ters are listed in Table |4] and their variability with time 
and X-ray flux is shown in Figure [8l These fits allow 
us to estimate both the maximum color temperature of 
the accretion disk and the total luminosity emitted by 
the disk and to compare these to IMBH and supercrit- 
ically accreting SMBH models in § 15.21 In this model, 
the accretion-disk emission contributes 30-60% of the to- 

Using XMM-Newton observation 0205590401 instead of 
0205590301 to obtain the power-law photon index and normaliza- 
tion does not affect the results. 



tal flux in the 0.3-10 keV band and the maximum color 
t emperature incr e ases w ith increasing flux. 

Stobbar t et all (|2006f ). who studied the highest S/N 
XMM-Newton spectra of 13 ULXs, concluded that nei- 
ther simple power law models nor simple MCD models 
provide adequate fits to the highest S /N spectra of ULXs. 
Better representations of their data are found using cool 
(~0.2 keV) or warm (^^1.7 keV) disk plus power-law mod- 
els, combination of cool black body plus warm MCD 
("dual thermal") models {bbody+diskpn in XSPEC), or 
MCD plus comptonized-corona models (diskpn+ eqpair in 
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Fig. 6.— NGC1365X1: MCD+power law fits to the 0.3-10 keV Chandra and XMM-Newton spectra. The MJD increases from top to 
bottom and left to right. The ObsIDs are given in the top right corner of each panel. The apparent residuals around 5-6 keV in XMM- 
Newton observation 0205590301 are likely an artifact of the background subtraction, since in this case the obscured active nucleus was 
particularly bright relative to XI in the hard band and this "feature" is absent from spectra extracted from regions with smaller extraction 
radii. 

2 keV. We therefore caution that the interpretation of the 
MCD plus power-law model fits of XI presented in § 15.21 
depends on the assumption that the hard band spectrum 
is well represented by a simple power-law arising in an 
optically thin corona, while the soft emission arises di- 
rectly from the accretion disk. 



XSPEC). The authors caution that depending on the 
choice of model, the inferred disk temperatures of their 
sample would cluster either at the "cool" or "warm" end 
of a luminosity-temperature diagram, and be interpreted 
as evidence of the presence of IMBHs (if cool) or SMBHs 
(warm). Even the best spectra (e.g., over 10'' counts with 
XMM-Newton) were inadequate to distinguish between 
these models. Thus, disk-temperature measurements 
based on spectra with less than a few thousand counts 
could give misleading results, especially when using the 
"wron g" spectral model. In addition, IStobbart et al.l 
()2006[) suggest that the evidence for a distinct curvature 
in the 2-10 keV spectra of ULXs (seen as a break in the 
power-law emission in 8/13 of their sources) can be used 
as a red flag against the interpretation of the hard band 
spectrum as emission from an optically thin corona. A 
broken power law model above 2 keV is not required by 
any of our individual observations nor by the joint fit to 
the highest S/N XMM-Newton and/or Chandra spectra 
of XI. Only the two longest XMM-Newton observations 
of XI, however, have >2000-counts spectra, which could 
be responsible for the non-detection of curvature above 



3.3. XI: Optical/IR Observations 

XI, whose X-ray position is known to better than 1" 
based on the Chandra images, was not detected on any 
of the optical or IR images available. Nevertheless, we 
use these images to obtain sensitive limits to the R, /, 
J, H, and K Johnson-Cousins band magnitudes of the 
ULX counterpart. The details of the lower- limit esti- 
mates are presented below and summarized in Table [5j 
We correct all optical and IR i mages astrometric ally us- 
ing the USNO-B catalog stars (jMonet et al.ll2003D in the 
immediate vicinity of the ULX source position. 

The left panel of Figure [i shows a 30 s VLT F0RS2 
(Very Large Telescope FOcal Reducer and low dispersion 
Spectrograph) R-Spccial filter^ observation of NGC 1365 

* The R-special filter is a Johnson-Cousins filter, slightly short- 
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Fig. 7. — NGC 1365X1; MCD+power-law fits to the 0.3-10 keV April 2006 Chandra spectra. The MJD increases from top to bottom 
and left to right. The ObsIDs are given in the top right corner of each panel. 



"T — I — I — I — I — I — I — r 








8 


o 

\ 





6 


■a 
o 





4 


n 
Ll_ 





2 










_l I I I I I I L_ 





> 



1.5 r 
1 - 
0.5 - 



52500 53000 53800 
MJD 




53850 



> 

^ 1.5 



O 

(U 
CL 



1 

0.5 







1 1 1 1 1 1 


r 1 




^ J 


-i- 






1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 ~ 



10 20 30 



40 



[10-1'* erg s-^ cm-2] 



Fig. 8. NGC 1365X1: variability of the spectral fit parameters for the MCD+power-law fits with MJD (left panels) and 2-10 keV band 
flux (right panels). 
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TABLE 5 

NGC 1365X1: Optical and IR Limits 



Telescope 


Filter 


T 


Aperture 


MJD 


ULX limit 


starE 


starW 


(1) 


(2) 


(sf 


(4) 


(5) 


(6) 


(7) 


(8) 


VLTFORS-2 


ij-special 


30 


1.5" 


52736 


R > 23.3 


18.5 ±0.3 


19.4 ±0.3 


HST WFPC2 


F814W 


10 


0.8" 


49732 


/ > 21.4 


16.81 


17.56 


HSTWFPC2 


F814W 


10 


0.8" 


49732 


/ > 21.6 


16.96 ±0.03 


17.71 ±0.05 


HST WFPC2 


F555W 


100 


0.8" 


51220 


V > 24.1 


19.8 ±0.1 


20.4 ±0.1 


2MASS 


J 




3.0" 


51602 


J > 17.2 


15.18 ±0.07 


16.32 ±0.17 


2MASS 


H 




3.0" 


51602 


H > 16.2 


14.56 ±0.06 


15.79 ±0.21 


2MASS 


K 




3.0" 


51602 


K > 15.9 


14.50 ±0.10 


15.32 ±0.20 



Note. — (1) and (2): Telescopc/instumcnt configuration and filter used in each observation; 
(3) Exposure time, T^xp, in seconds; (4) Aperture radius used to measure source fluxes; the 
background was measured locally within an aperture with twice the area of the source aperture; 
(5) Modified Julian Date (MJD) of the observation; (6) ULX limit in the standard Johnson- 
Cousins filters; (7) and (8) Magnitudes of the two stars (to the East and West of the ULX 
position, respectively; see Figure [9j used for magnitude calibration in each Johnson-Cousins 
filter. 




Fig. 9.— The optical region around NGC 1365X1 through the VLT-FORS2 R-band filter (left) and the HST F555 (top right) and F814 
(bottom right) filter s. T he ULX is marked by 2"-radius circles on all images; 1.5" (0.8") circles denote the W and E stars on the VLT 
(HST) images (see 13.31 1. The intensity of the VLT image is scaled to reveal the diffuse extended emission at the position of the ULX as 
well as the two faint point sources 3.7"S and 4"NE of the X-ray position of XI. 
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within 50" of the X-ray position of XI. SmaUer region 
cutouts around the XI X-ray position in two of the HST 
filters are shown in the right panel of Figure ^ The 
ULX X-ray position coincides with a faint spiral exten- 
sion; two faint sources are seen 3.7"S and 4"NE of the X- 
ray position, respectively. Considering the large separa- 
tions, neither of these sources is a likely optical counter- 
part, but the optical counterpart of the ULX was clearly 
fainter than both sources during the observation. The 
two bright stars - 16"E and 13"W of the ULX in Fig- 
ure [9]have USNO-B 1.0 R-band magnitudes of 18.5 ± 0.3 
and 19.4 ± 0.3, respectively. Judging by the flux ratios 
between those stars and the faint sources ~ 4" S and 
NE from the ULX X-ray position, the faint sources have 
R-band magnitudes of ~ 22.4 and ^ 22.3, respectively. 
We chose a local background region a few arcsec NW 
of NGC 1365X1 (but still within the diffuse spiral arm 
emission) and aperture radii of 1.5" for all flux measure- 
ments. The R-band magnitude limit for the ULX emis- 
sion is i? > 23.3. 

Two HST WFPC2 observations of the central region of 
NGC 1365 (see the right panel of Figure[9]) were taken in 
1995 with the F814W (IDs) and F555W filters (100s). 
The HST F814W filter is close to the Johnson-Cousins 
J-band. Using the /-magnitudes of the same two bright 
stars (E and W of the ULX source position in Figure [9]) 
as photometric calibrators, we estimate a lower limit to 
the apparent magnitude at the ULX position of / > 21.4 
(see the 2d row of Table [5]). If we use the indepen- 
dent HSTWFPC2 photometric calibration,^ we obtain 
the values given in the 3d row of Table O and a ULX 
source position lower limit of / > 21.6, in agreement 
with the lower apparent-magnitude limit obtained using 
bright-star calibration. The lower limit to the V-band 
apparent magnitude at the ULX position given in the 
4th row of Table [5] was estimated in similar fashion using 
the 100 s F555W-filter HST exposure. 

2MASS coverage of the region around the ULX source 
position is also available. We use relative photometry 
(using the same bright stars E and W of the source) 
to estimate the H, J, and K band lower limits. The 
results, using 3"-radius source aperture measurements, 
appear in the last three rows of Table [5l We comment 
on the implications of all optical/IR apparent magnitude 
limit results in § 15.1.21 

4. NGC1365X2 

NGC 1365X2 (03:33:41.85, -36:07:31.4) is a highly vari- 
able X-ray point source in the north tip of the eastern 
spiral arm of NGC 1365 (see Figure [2]). Based on the 
highest luminosity reached during one of the April 2006 
Chandra observations, which exceeds 4 x 10'*°ergs~^ in 
the 0.5-10 keV band, the object is among the most lumi- 
nous ULXs known. 

We extracted spectra for the X2 source from all Chan- 
dra observations as described in § 13.21 for NGC 1365X1. 
Because of the proximity of X2 to other X-ray point 
sources (the closest one is only 9"S of X2, the second 
closest 20"W) combined with the larger XMM-Newton 
point spread function and the faintness of X2 during 
the XMM-Newton observations, we do not attempt to 
extract spectra/upper limits from these exposures. We 

^ http:/ /www. stsci.edu/documcnts/dhb/wcb/c03_stsdas. fm3.html 
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Fig. 10. — NGC 1365X2: X-ray spectroscopic fits to the high- 
est S/N observation, Chandra Obs.ID. 6870. All fits include an 
intrinsic absorber withJVn.i = 6 — 9 X 10^^ cm~^ as well as Galac- 
tic absorption (see i) 15.2.21 for details). Top: Power-law fit (inset) 
and 3(T confidence contours for the power-law photon index and 
the column density of the intrinsic absorber. Middle: MCD fit 
(inset) and Sa confidence contours for the normalization and the 
maximum color temperature. Bottom: The kerrbb fit (inset) and 
3(T confidence contours for the black hole mass and accretion rate. 
The spectral fits, shown as insets, give the data overlaid with the 
model in the upper panels and the fit residuals in the lower panels. 

attempted to fit the Chandra spectra with both power- 
law and simple MCD models. Assuming galactic absorp- 
tion, we find the simple power-law or MCD fits are un- 
acceptable and intrinsic absorption is necessary in both 
cases. The absorbed power-law fits are summarized in 
Table [6l Due to the small number of counts, we use the 
Cash-statistic fit for all but one of the spectral fits - the 
highest signal-to- noise (S/N) spectrum (Obs.ID. 6870), 
which has 535 counts in the 0.3-10 keV band, was binned 
to 15 counts per bin and fit using statistic. We show 
the power-law fit and the 3a confidence contours for the 
photon index and the column density of the intrinsic ab- 
sorber for the highest S/N fit of X2 in the top panel of 
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Fig. 11.— The optical region around NGC 1365X2 through the 
VLT-FORS2 R-band filter. The X-ray point sources are denoted 
by 3"-radius circles. X2 coincides with the extended emission of 
an HII region. 

Figure [TOl For all spectra with less than ~ 200 counts 
in the 0.3-10 keV band, we fix the power-law spectral in- 
dex to r = 1.3, the value obtained for the two spectra 
with the largest numbers of counts. A cold absorber of 
iVn i = 6 — 9 X 10^^ cm^^ is required by all fits, including 
the MCD fits. 

The absorbed MCD fits are statistically indistinguish- 
able from the absorbed power-law fits presented above. 
For all April 2006 Chandra observations the maximum 
color temperature estimated from the MCD fits was 
around 2keV, with large uncertainties (Itr uncertain- 
ties ~ 0.3 — 0.6 keV). The 3a confidence contours for 
two interesting parameters for the highest S/N MCD fit 
are shown in the middle panel of Figure [TOl The ad- 
dition of a comptonized component (e.g. a power law) 
does not improve the fits. We caution that, consider- 
ing the small number of total counts and even smaller 
number of counts above 5 keV, it is conceivable that the 
absorbed MCD model is inappropriate, and the disk- 
temperature constraints are unreliable. ULXs and X-ray 
binaries whose MCD fits result in unusually high disk 
temperature s similar to those obtaine d for X2 have been 
known since iMakishima et al.l (|2000( l . iMakishima et al.] 
(|200Cll ) suggested that the unusually high temperatures 
are naturally explained by emission from a disk around 
a rotating black hole. The bottom panel of Figure [10] 
shows a multi-temperature blackbody model fit to the 
highest S/N X2 spectrum which assumes a thin, steady- 
state, general relativist ic accretion di sk around a rotat- 
ing (Kerr) black hole (jLi et all [20051 ). Both the MCD 
and the kerrbb fits shown in Figure [TOj are acceptable 
ix^/DoF = 25/31 and x^/DoF = 24/32 , respectively) 
and we comment on the results in § 15.2.21 

X2 f alls within a 1995 ROSAT HRI observation. ISinghl 
(|1999D claim a 3.5cr detection with 0.0020 ± 0.0006 cps 
(■Po.5-iokev ~ 1-3 X 10~^^ ergs^^ cm^^ assuming a power 
law model with F = 1.3 and TVh.i = 9 x lO^^ cm-^). We 
reanalyzed the ROSAT HRI observation and find only a 
very faint source with 0.00053 ± 0.00026 cps at that po- 
sition (~ 2cr); we use this as an upper limit to the X-ray 
count rate of X2 at this epoch, arriving at an obscuration 
corrected fiux of Fo.5_iokGV < 3.5xl0^^'^ergs~^ cm^^ for a 
power law model with F = 1.3 and N^a — 9 x 10^^ cm~-^. 



X2 is one of the three sources that were not detected dur- 
ing the 2002 Chandra observation with an upper limit of 
CR< 0.0004 cps or Fo.5-iokov < 0.9 x lO-^ergs"! cm'^. 

5. IMBHS, SUPER-EDDINGTON HMXBS, OR 
BACKGROUND AGNS? 

Super- or hyper-novae can reach X-ray luminosities of 
over 10^°ergs~^, but they are expected to be steady 
sources on short timescales and to decline in luminos- 
ity slowly over periods of months to years. The strong 
variability of both NGC 1365X1 and NGC 1365X2 over 
extended periods of time (factors of 10 over years) and 
their short time flares argue against supernova emission 
and imply a compact accreting source in both cases. In 
addition, neither of the two ULXs discussed here has the 
thermal emission-line dominated spectrum from shock- 
heated plasmas which is characteristic of supernovae. 
Below we discuss the IMBH, super-Eddington SMBH, 
and background AGN possibilities for the two NGC 1365 
ULXs, focusing on their spectral characteristics, lumi- 
nosity variability, and optical/IR-to- X-ray ratios. 

5.1. Background AGNs? 

We argue that both NGC 1365X1 and NGC 1365X2 
are unlikely to be background AGNs. Based on their po- 
sitions - XI coincides with one of the diffuse weak spiral 
arms south of the NGC 1365 nucleus, while X2 coincides 
with an HII region at the tip of the main NE spiral arm, 
as well as the negligible number of background sources as 
bright as these ULXs expected for a similar-sized random 
patch of sky, XI 's unusual X-ray-to-optical/IR ratios and 
its small intrinsic X-ray absorbing column, an IMBH or 
a supercritically accreting SMBH are more likely in both 
cases. We cannot rule out that XI, which has no optical 
counterpart, is not an optically/IR obscured AGN with 
little X-ray obscuration with certainty, but the associa- 
tion of NGC 1365X2 with the H II region in the NE arm 
of the galaxy makes it extremely unlikely that the source 
is a background AGN. 

5.1.1. log N — log S Argument 

According to th e log TV — log S Lockmann Hole results 
of iHasinger et al.l ()2001f ). we expect < 0.2 sources per 
square degree with Fo,5_iokcv > 5 — 6 x 10"^'^ ergss"^. 
Consequently, < 2 x 10~^ sources with flux similar or 
higher than that of the maximum flux of XI or X2 are 
expected within 2' of the center of NGC 1365. This, to- 
gether with the fact that XI is located in one of the faint 
spiral-arm extensions of NGC 1365 and X2 coincides with 
an H II region at the tip of the NE arm, make it unlikely 
that either of the sources is a background AGN. 

5.1.2. X-ray-to-Optical/IR Ratios 
NGC 1365 XI 

As noted in § 13.31 a-ny optical counterpart of the XI 
source is fainter than R > 23.3. This limit is con- 
sistent with any stellar mass companion, including the 
brightest supergiants. Computing the X-ray-to-optical 
ratio, log[X/0] = logJ^x + OAR + 5.61, we find lower 
limits between 1.5 and 2.7 for the full range of X-ray 
fluxes observed for XI. Typical X-ray-to-optical ratios 
for AGNs range between —1 < log[X/0| < 1 (e.g., 
iMaccacaro et al.l 119881 : iHornschemeier et al.l l2001f) , with 
the upper range corresponding to blazars. About 90% 
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of all AGNs have log[X/0] < 2. Since the observed 
lower limit to the X-ray-to-optical ratio for XI is smaller 
than this value for only 3 of the 12 X-ray observations, 
we consider it highly unlikely that this ULX is a back- 
ground AGN. For comparison, the 8 ULXs which were 
found to be b ackground AGNs (or a narrow l i ne gal axy in 
one case) by [Gutierrez fc Lopez-Corredoiral ()2007D have 
0.04 < log [X/0] < 0.6. The two additional X-ray point 
sources from Table [1] with optical counterparts in the 

USNO-Bl.O catalog have log [X/0] 0.4 (X17) and 

log [X/0] ~ 0.3 (X15). X15, whose optical counterpart 
is a point-source, is likely a background AGN. 

We estimated a limit of / > 21.6 (F/ < 6^Jy at 
0.9 Aim) for the /-band magnitude of NGC 1365X1. The 
X-ray-to-I-band ratio, log[X//] = log [Fq 5_2kev/^/], has 
a lower limit ranging between —0.1 and 1.7 for the 
full range of observed 0.5-2 keV fluxes, with a most 
likely lower limit between 0.4 and 0.7, indicating that 
the source could be a very rare optically faint/X-ray 
bright background A GN at a redshift >2 (see Figure 4 of 
iBrandt fc Hasinger, 20051. If that was the case, however, 
the measured flux would imply an X-ray luminosity of 
~ lO'^^ ergs"^. According to the X-ray luminosity func- 
tion of lHasinger et al.) (|2005f ) type 1 (unobscured) AGNs 
with such high luminosities are rare: ~ 10~^ Mpc~'^. The 
analysis of the X-ray background (e.g.. lGilli et al.ir2007f l 
further shows that the ratio of obscured to unobscured 
AGNs is < 1 for luminous AGNs, implying that there are 
< 200 AGNs in total of comparable luminosities between 
redshifts of 1.9 and 2.1 in the whole universe, which to- 
gether with the large variability of XI, uncharacteristic 
for luminous AGNs, make it very unlikely that XI is a 
luminous redshift >2 AGN. 

Using the 2MASS observations, the Ks > 15.9 limit of 
XI corresponds to Fk^ < 0.29 mJy at 2.2 /im, indicating 
a lower limit of the X-ray-to- i^Tg flux ratio of ^0.2 (rang- 
ing between .04 an d 0.6 for the observed X-ray range). 
IWilkes et al.l (|2002D studied the X-ray properties of red 
2MASS AGNs. According to their Figure 1, typical red 
2MASS AGNs have X-ray-to- i^g flux ratios < 0.05, while 
the z <1 broad- line AGNs form lElvis et all (p^994) have 
X-ray-to- ifs hux ratios between 0.05 and 0.7. The lower 
limits to the X-ray-to- /sTg flux ratio for NGC 1365X1 is 
consistent with an unobscured broad-line AGN, but this 
interpretation is inconsistent with the observed optical 
limits. 

NGC 1365X2 

In the optical image, X2 is associated with a bright 
extended emission region in the eastern spiral arm of 
NGC 1365 (see Figurefll]), 86" S from the NGC 1365 nu- 
cleus (equivalent to ~ 8 . 8 kpc) , which was classified as 
an H II region by 'Hodgd (|1969h . Using the USNO-B 1.0 
catalog's i?-band magnitudes of the optical counterpart, 
Rl = 13.25 and R2 — 13.75, we estimate an X-ray-to- 
optical ratio of —2 < log [X/0] < —1, depending on the 
X-ray flux level. The H II regio n has been studied in de- 
tail bylRov fc Walsh ' (1997^ and'Ravindra nath fc Prabhul 
((200l . It is a typical large star-forming region with a 
diameter of 670pc, Mv = —12.4, and an Ha line lu- 
minosity of 1.6 X lO'^^ergs"^ (S. Rav indranath, private 
communication). iRov fc WalshI ()1997f ) study the optical 
spectrum of the H II region (number 33 in their Table 2) 
in comparison with other H II regions in NGC 1365. The 
observed narrow-line ratios suggest an ionization-bound 



nebula which is photoionized by massive stars. There 
is no evidence that the X-ray emission of NGC 1365X2 
affects the observed line ratios, which is consistent with 
its interpretation as an X-ray binary (whose energy in- 
put would be too small to affect the observed narrow 
line-ratios of a giant nebula). 

According to the X-ray observations, X2 contributes 
about 5 X lO'*^ ionizing photons per second, and the 
Ha luminosity emitted as a consequence of this ioniz- 
in g flux should be about 5 x lO^^ergs"^ (see Section 5.3 
of iSoria et al.l [2005, and references therein for the de- 
tails of the computation). This is a factor of 30 fainter 
than the measured Ha line luminosity of the H II region, 
confirming the conclusion that X2 does not contribute 
substantially to the photoionization of the surrounding 
nebula. 

5.2. L-T Diagram Results 
5.2.1. NGC1365X1 

Panel A of Figure [12] shows the Luminosity- 
Temperature {L-T) diagram for NGC 1365X1. The 
estimated MCD temperatures and accretion disk- 
luminosities follow a, L (X T" relation, with a — 
Panel B of Figure [T^] displays the theoretical tracks 
for 20, 40 and 60 M© black hole primaries, accret- 
ing at sub-Eddingt on (colored red) or sup er-Eddington 
rates (colored blue; iPoutanen et al.l [20071) . The XI re- 
lation agrees with the L oc relation expected for 
a black hole primary with mass between 40 Mq and 
60 Mq, but the observed maximum luminosities during 
the Chandra flare of Aprfl 2006 {L^ > 10'^° ergs' ^) 
require super- Eddington accre ti on ra tes and/or beam- 
ing. In the Pou tanen et al.l (|2007l) model, the ob- 
served luminosities can reach L « LBdd(l + 0.61nrh) 
where ih is the accretion rate in units of Edding- 
ton.^" The maximum color temperature of the disk is 
= x 1.6 X (m-i/4) X (1 - 0.2m-i/3)keV, where 
m is the black hole mass in units of Mq, and f^ is 
the color-correction factor, — 1.7 (r„,ax is the max- 
imum temperature o f the advective disk, see eqn. 36 of 
IPoutanen et al.ll2007f ). A 40-60 Mq black hole primary is 
consistent with all XMM-Newton and Chandra points in 
Figure lT^ within 3cr of their estimated uncertainties if the 
accretion rate varies between .1 and 16 of the Edding- 
ton rate in the IPoutanen et al] (^007) model. The mea- 
sured temperatures during the Chandra flare, are sys- 
tematically high relative to the super-Eddington 60 Mq 
model (by about 30% or ^ 3ct). A variation of the color- 
correction factor with accretion rate for l<rh<16 could 
account for this discrepancy. 

We note that a less massive primary is problematic not 
only because it fails to explain the observed temperature- 
luminosity pairs in Figure I12[ which could, as noted in 
§ 13. 2i be affected by the specific model fit assumptions. 
A 10 Mq black hole primary, for example, is also prob- 
lematic because the required beaming factors or super- 
Eddington accretion rates (a beaming of over a factor of 
30 or, alternatively, > 10~^ Mq available for accretion in 
^2 days during the Chandra flare in the Poutanen et al. 
2007 model) are uncomfortably large. 

In the supercritical accretion model of lOhsuga et all (120051 ), 
the luminosity can exceed the Eddngton luminosity by a factor of 
up to ~10 for face on view and rh=100-1000. 
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Fig. 12. — NGC1365X1 - the multi-color disk (MOD) temperature vs. the accretion disk luminosity: data {Panel A) and theoretical 
models {Panel B). The XMM-Newton/ Chandra points are denoted by open/solid circles. A power-law fit to the estimated MOD parameters, 
L oc T°' , with a = 2.8^Q g, is given with the dashed lines in both panels. The blue lines show the (color corrected) maximum temperature 
(Tlnax) of supercritically accreting black holes with M=2OM0, M=4OM0, and M=60Mq l lPoutanen et al.ll2007l '). The accretion rate varies 
between 1 < rh < 800 for T„ax (each triangle corresponds to a doubling of the accretion rate). The red lines show the L oc relations 
expected for 20 Mq, 40 Mq, and 60 Mq black holes with accretion rates 0.1 < rfa < 1 (each dot corresponds to an increase in the accretion 
rate by 0.1). 



We conclude that the L-T diagram of XI is consistent 
with the expectations for a 40-60 M© black hole primary 
accreting at 0.1<m<16 from a massive companion star, 
offering no support for the IMBH hypothesis. 

5.2.2. NGC 1365X2 

As described in § 21 the X2 spectra are consistent with 
absorbed power-law, or absorbed MCD fits, and the cur- 
rent data is of insufficient quality to distinguish between 
these or more complex models (e.g. models including a 
combination of a MCD and a comptonized component) . 
The MCD fits to the X2 spectra are inadequate to con- 
strain the the maximum color temperature and inner ra- 
dius of emission for each individual observation and as- 
suming a simple MCD model could lead to unreliable 
temperature estimates. As a result, the following inter- 
pretation of the results depends critically on the assump- 
tion that X2 is in a disk-dominated state and the MCD 
or kerrbb models are appropriate. 

Assuming that X2 is in a disk dominated state and 
the MCD model is appropriate, the disk temperature 
estimates are too high for a sub-Eddington accretion 
around a non-rotating black hole (r^,„,ax < l.TkeV). The 
iPoutanen et all ()2007l ) model with a SMBH primary can 
reproduce the observed disk temperatures for M < 5 Mq , 
but cannot reproduce the maximum luminosity (a fac- 
tor of >50 higher), since in that model the luminosity 
exceeds Eddington by a factor of 1 -I- 0.61nrii, which 
is less than a factor of 5 for rh < 10'^. Other su- 
percritical accretion models (e.g. lOhsuga et al.l l2005f ) 
also predict luminosities which exceed the Eddington 
by factors of a few, up to 10-15. If we would believe 
the high accretion disk temperatures, iMakishim a et al.] 
jlpOO) have suggested a solution - a rotating black hole 
would naturally explain the higher disk temperatures. 
For X2 a maximally rotating black hole (a = 0.998) 
with disk inclination fixed to 6* = 85 (close to edge 
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Fig. 13.— NGC 1365X2 luminosity variation with time. The 
Chandra data is shown with solid circles, the _RO>Sj4T-HRI limit 
with an open circle; the arrows indicate upper limits. 



on), we obtain M — 210 ± 70 Mq with an accretion rate 
M = (1.8 ±0.3) X 10^0 gs^i (rii w 0.4 - 0.5), and an in- 
trinsic absorption of A^H,i = 6±2 x 10^^ cm~^. The bolo- 
metric luminosity of this model is ~ 1 — 2 x 10""^ ergs^^, 
similar to the kerrbb-&t value of 3.5^^ x lO'*" ergs~^. 
Consequently the X2 spectra are consistent with emis- 
sion from a disk around a maximally rotating black hole, 
for a large range of IMBH masses, M 80-500Mo. The 
S/N however is not sufficient to distinguish those mod- 
els from simple absorbed power-laws, and the evidence 
provided by this spectroscopic fit is only suggestive. 

5.3. Implications of the X-ray Variability 
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TABLE 6 

NGC 1365X2: Absorbed Power-Law Fits 



ObsID 


-^2-10 koV -^0.5-2 kcV 


r 


Nh 


Data/Method 


stat/DoF 


(1) 


(2) (3) 


(4) 


(5) 


(6) 


(7) 



3554 


< 0.7 


< 0.2 


1.3 fixed 


9 fixed 


ACIS-S 




6871 


5 fi+3-8 

°-"-3.3 


1 7+0.6 
1- '-0.6 


1.3 fixed 


7±3 


ACIS-S/C-stat 


40/44 


6872 


42+^ 


12+" 


1.3 ±0.2 


9±2 


ACIS-S/C-stat 


165/228 


6873 


37+1 


-,,+2.0 
1^^-1.6 


1.3 fixed 


10 ± 1 


ACIS-S/C-stat 


175/187 


6868 


25+1 


7 4+1.9 


1.3 fixed 


9 ± 2 


ACIS-S/C-stat 


122/151 


6869 


12tl 


3 6+^ '^ 

J.D_j 2 


1.3 fixed 


8±1 


ACIS-S/C-stat 


73/102 


6870 


63+_\l 


20+'''* 
^'-'-12.4 


1.3 ±0.2 


9 ± 2 


ACIS-S/x^ 


26/31 



Note. — The fiuxcs arc quoted in units of lO^^" crgs^^ cm^^. (1) Observation 
ID; (2) Absorption-corrected fiux in the 2-10 keV band; (3) Absorption-corrected 
fiux in the 0.5-2 keV band; (4) Power-law photon index; (5) Intrinsic absorption 
column density in units of 10^" cm~^ ; (6) Data and fitting method used; (7) / DoF 
or C-stat/DoF of the fit. 



The X-ray spectra of classical X-ray binaries initially 
get softer as their flux increases, which is not observed 
for either NGC 1365X1 or X2. If anything, the X-ray 
spectrum of NGC 1365X1 gets harder in the higher flux 
states. Other ULXs are known to show this behavior on 
shorter timescales which was suggested to arise in cool- 
accretion disk plus optically thick corona models of ULX 
emission as the coro na is heated in the h igh flux states 
(e.g. NGC 5204X1; [Roberts et aP [2006h . iSoria et all 
()2007t l interpret the hardening with increased flux as the 
effect of an additional ionized-absorption component - a 
flux decrease in the 0.5-1.5 keV region which looks like a 
continuum hardening in power-law fits with neutral ab- 
sorption. 

From Figure m NGC 1365X1 appears to be persis- 
tently luminous in the X-ray band (typically a few 
times lO^^ergs"^) over decades with occasional short 
(~days), but not very powerful (< a factor of 10) flares. 
NGC 1365X2, on the other hand, has more transient be- 
havior - it was not detected before April 2006, when it 
suddenly reached unobscured luminosities of a few times 
]^q40 gj.gg-1^ varying on a timescale of days (Figure [T5)) . 
Assuming either a SMBH or an IMBH binary, the per- 
sistent high state of XI can be interpreted as evidence 
of a massive donor star (>2O-3OM0) supplying a large, 
stationary, but not necessarily steady (considering the 
Chandra flare on timescale of days) accretion disk. As 
noted in ij 15.21 above, the 2006 Chandra flare of XI is too 
small to be associated with a disk instability, as observed 
in SXTs in LMXBs, which result in accretion rate and 
l uminosity increa s es of factors of ^ lO''. 

iKalogera et al.l (2004) studied the long term tran- 
sient behavior of SMBHs and IMBHs due to thermal- 
viscous disk instability for ULXs in young stellar pop- 
ulations. For the most likely donors in such environ- 
ments, with masses in excess of 5Mq, they inferred a 
minimum black hole mass for transient behavior, which 
is > 50 M0. Based on this argument and the observed 
variability patterns of the two ULXs, it is more likely that 
NGC 1365 X2, which is a transient source associated with 
an H II region, contains an IMBH. 

5.4. Free-Floating Supermassive Black Hole? 

Recent numerical relativity simulations of back hole 
mergers suggest that large gravitational recoil veloci- 



ties of up to several 1000 km s~^ are possible for merg- 
ers of maximally spinning, equal-mass, bl ack holes with 
anti- al igned spins in th e orbit al plane pCoppitz et al.l 
20071: iCampanelli et all 120071: [Gonzales et all 120071: 



Herrmann et al.l 120071 : iBaker et al.ll2007t) . Alternativelv. 



a gravitational slingshot of three or more supermassive 
black holes can re sult in the ejection of one of them 
from the nucleus ()Hut fc ReesI 1199^ iHoffman fc Loebl 
|2007() . Luminous off- nuclear X-ray sources are poten- 
tial counterparts to the s e recoiling black hole s (e.g., 
iMadau fc Quataert) [200l iLibeskind et al.l I2OO60 . Is it 
possible that either XI or X2 is a "free-floating" super- 
massive black hole which was expelled from the center 
following a merger with a large recoil velocity or a grav- 
itational slingshot in the past? 

A recoiling BH with a kick velocity below the es- 
cape velocity will oscillate within the galaxy with a 
spatial amplitude that depends on the kick v elocity 
(jiviadau fc Quataertj 12004 iGualandris fc MerrittI [2008) . 
Oscillations typically damp within a time scale of ^ 
10^^^ years, which is comparable to the time scale of 
a galaxy merger. At present, we do not see evidence for 
the occurence of a recent merger in NGC 1365. In partic- 
ular, a nearly equal-mass merger would have resulted in 
the formation of a large bulge and randomized the central 
velocities, contrary to the observed strong bar pattern in 
NGC 1365. 

Furthermore, within the recoil scenario, we would not 
expect a central black hole in NGC 1365 at the cur- 
rent epoch. However, NGC 1365 does harbour an ac- 
creting supermassive black hole at its core, as is evi- 
denced by the presence of broad Balmcr li nes and A GN- 
typical emission-line ratios (e.g.. .Schulz et al.lll999D and 
by th e detection of rapid X-ray variabilitv ([Risaliti et al.l 
|2007[) . In addition, the coincidence of X2 with an H II re- 
gion would be hard to explain under the recoil/slingshot- 
ejection scenarios. With respect to XI and any other 
off-nuclear point source not associated with a dense re- 
gion, assuming that its black hole mass does not ex- 
ceed 10^ Mq, additional evidence against the free-floating 
black-hole interpretation comes from the availability of 
accretion material. To reproduce the observed luminosi- 
ties with a 10^ Mq black hole (assuming a temperature 
of 10"* K) we would need to embed it in a medium with 
density of at least n « lcm~'^, which is typical for a 
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normal HII region, not the diffuse spiral arm where XI 
is located. For the typical interstellar medium which has 
higher temperatures and lower densities (n ~ 10"'^ cm~'^ 
and T ^ 10^ K), the Bondi accretion rate is a factor 
of 10^ lower. While ejected black holes can be tem- 
porarily fueled by stellar tidal disruptions or stellar mass 
loss from the stars which a re bound to them upon recoil 
(jKomossa k. Merrit3l2008[ ) , these events are rare and are 
very unlikely to operate in two sources simultaneously in 
one nearby galaxy. Based on all of these arguments, none 
of the off-nuclear X-ray sources of NGC 1365 appears to 
be a compelling recoil/slingshot-ejection candidate. It 
is still possible, however, that recoils and/or black holes 
ejected by a gravitational slingshot hide among the ULX 
population as a whole. 

6. SUMMARY AND CONCLUSIONS 

We present 26 X-ray point sources detected in a se- 
ries of Chandra observations of the giant spiral galaxy 
NGC 1365. The cumulative distribution of sources 
brighter than a given X-ray flux is compatible with pre- 
dictions based on the star-formation rate of the galaxy, 
as expected if high-mass X-ray binaries dominate the 
extended X-ray emission of star-formin g galaxies (e.g., 
iGilfanov et an[200l : iGrimm et al1[200l . 

Two of the X-ray point sources, NGC 1365X1 and 
NGC 1365X2, are highly ultraluminous, with peak lu- 
minosities exceeding 3-5xl0'*''ergs~^. Both ULXs are 
most likely accreting compact objects (with SMBH or 
IMBH primaries), since both the super-/hyper- nova ori- 
gin and the background or "free floating" AGN interpre- 
tations appear highly unlikely. Both sources are among 
the most luminous ULXs known and as such, among the 
most likely candidates for IMBH binaries. Their lumi- 
nosities and optical counterparts (a lower limit for XI 
and a giant H II region for X2) are consistent with an ac- 
creting IMBH in both cases; assuming no beaming and 
accretion rate equal to the Eddington rate (computed 
from the highest Chandra luminosities), the black hole 
masses are ~ 250 in the case of XI and ~ 340 M© in 
t he case of X2. 

iRemillard fc McClintocfl ()2006f ) define three outburst 
states of SMBHs, based on data of galactic BH bina- 
ries: (1) the thermal or accretion-disk dominated state 
(with disk fraction of > 80%; previously known as the 
High/Soft state), (2) the hard (1.4 < T < 2.1) or power- 
law dominated state (the Low/Hard state), and (3) the 
steep power law (SPL) state, in which both disk emission 
(typically <50%, but up to 80%) and a steep power law 
(r > 2.4) contribute (the Very high state). Active nuclei 
with supermassive black holes are typically observed in 
the hard state, in which case a radio jet is often present 
in the case of X-ray binaries as well as low-luminosity 
AGNs. It is not clear whether accreting IMBHs will ex- 
hibit the three states observed in galactic BH binaries. 

The observed X-ray spectra for XI and X2 are consis- 
tent with power-law models including intrinsic absorp- 
tion, as well as thermal accretion-disk (plus a power 
law in the case of XI, which contributes about 50% of 
the flux) emission models. The 12 Chandra and XMM- 
Newton spectra of XI are equally well fit by power-law 
emission modified by intrinsic neutral absorption and a 
multi-color disk (MCD) emission plus about equal con- 
tribution from a comptonizing corona (reminiscent of the 



SPL/ Very- high state). The MCD fits to the X-ray spec- 
tra of XI imply maximum color temperatures which are 
too high for the black hole masses inferred from the max- 
imum black-body luminosity {'^ 100 M©), if we assume 
isotropic emission and sub-Eddington accretion rates, 
suggesting that beaming and super-critical accretion also 
play a role. Assuming that the MCD plus power-law 
fits provide an accurate representation of the spectra, a 
black- hole mass primary in the 40-60 M© range is con- 
sistent with the observed L-T diagram for XI. 

In the case of X2, the observed X-ray spectra are 
consistent with intrinsically absorbed power-law, MCD, 
and maximally-rotating IMBH fits and the insufficient 
photon statistics prevents us from obtaining reliable 
accretion-disk temperature estimates. The MCD tem- 
perature estimates obtained (~2 keV) are unrealistically 
high for even the super-critically accreting SMBH mod- 
els, suggesting that the simple power-law model or the 
rotating black hole model provide more viable explana- 
tions for the observed X-ray spectra. A rapidly-rotating 
(a - 0.998) IMBH with M =80-500 Mq is consistent 
with the highest S /N Chandra spectrum of X2 and can 
reproduce the observed luminosities. An IMBH inter- 
pretation for X2 i s also consistent with the theoreti- 
cal expectations of iKalogera et all (|2004!), who inferred 
a minimum black hole mass for transient behavior of 
M > 50 M0 for a binary with a donor star more mas- 
sive than 5 M0 . 

Absorbed power-law models provide equally good fits 
to the X-ray spectra of XI and X2. The results of 
the power-law model fits for XI and X2 are reminis- 
cent of the hard (formally, Low/Hard) state of galac- 
tic BH binaries. The power-law slopes deduced for XI 
(1.6 < r < 2.2) and X2 (P = 1.3 ± 0.2) are con- 
sistent with those of galactic BH binaries in the hard 
state, where the comptonized component dominates; the 
observed luminosities are however very high, suggest- 
ing that super-Eddington accretion modes and or beam- 
ing/coUimation also play a role. The maximum lumi- 
nosities reached (3-5xl0^°ergs~^) are factors of 10-20 
higher than those of the < 20 M© black- hole X-ray bi- 
naries studied in our Galaxy, pushing the limits of the 
current supercritical accretion models and requiring very 
high beaming/coUimation. 

The most likely black-hole masses are probably around 
40-60 Mq in the case of XI, and 80-500 M© in the case 
of X2, although most of the arguments for such high 
masses come from the thermal components of the spec- 
tral fits, which could be unreliable. We consider masses 
of < 20 Mq less likely also because they are harder to 
reconcile with typical expectations for the supercriti- 
cal/beamed emission models available^ 

While this paper was in preparation. ISoria et al.l ()2007[ ) 
published an analysis of NGC 1365X1 and concluded 
that the spectra are best fit by absorbed power-law com- 
ponents in all cases, and an accretion disk component 
contributes at most 10-25% of the emission of XI in se- 
lect higher S/N spectra. The authors suggest that this 
is compatible with a « 50 — 150 Mq mass black hole, 
where the accretion disk is colder (0.2-0.4 keV) and the 
inner parts of the accretion flow are in the form of an 
outflow or a Comptonizing corona, which is responsible 
for the power-law emission. The current data, despite 
its high photon statistics, is insufficient to distinguish 
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between different scenarios without the assumption of a 
specific theoretical model (e.g., power- law emission with 
or without contribution from the accretion disk and in- 
trinsic absorption) followed by a targeted interpretation 
of the results. Higher S/N X-ray spectra and better con- 
straints on the short term variability of XI, as well as 
better theoretical understanding of the most luminous 
ULXs (which will allow us to select the proper theoreti- 
cal models to fit), will allow us to distinguish those sce- 
narios in the future. It is reassuring that despite our 
different interpretations of the X-ray spectra of XI, we 
reach similar conclusions about the possible black hole 
mass of XI. 

The optical counterpart of X2, a luminous HII region 
in the spiral arm of NGC 1365, is the most likely place to 
find an IMBH - both in terms of formation mechanisms 
for IMBHs, which require very dense environments, and 
in terms of the availability and capture of a suitable 
donor star to supply the accretion material. The lack of 
an optical counterpart of XI makes the association with 
an IMBH more problematic, but not implausible; recent 
theoretical work suggests that for primordial IMBHs (one 
of three possible formation mechanisms, all of which re- 
quire a dense environment), most IMBHs with masses 
< 1000 M0 would be ejecte d from the dense globu- 
lar cl uster where they formed llHolley-Bo ckclmann ct al.' 
|2007( ). Consequently, the lack of an optical counterpart 
for XI does not automatically preclude an IMBH pri- 
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